REPORT  DOCUMENTATION  PAGE 


Form  Approved  0MB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 
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4.  Statement  of  the  problem 


Flexoelectric  effect,  defined  as  a  linear  coupling  effect  between  a  strain  gradient  and  an  induced 
electric  polarization  [1],  has  drawn  scientific  interests  in  recent  years  [2]-[14].  This  phenomenon 
is  generally  expressed  in  tensor  form  as  [15], 


Direct :  P;  = 


ds. 


o 


dXi, 


(1) 


Converse  : 


ijkl 


dx, 


(2) 


where  P;,  and^Cnj  are  the  induced  polarization,  flexoelectric  coefficient,  elastic  strain, 

and  axis  of  the  gradient;  S^j,  and  are  the  strain,  converse  flexoelectric  coefficient 
associated  with  strain,  and  electric  field,  respectively. 


In  principle,  the  flexoelectric  effect  is  a  universal  phenomenon  exhibited  by  all  dielectrics  [16], 
Consequently,  non-piezoelectric  materials,  such  as  cubic  materials  and  isotropic  materials,  still 
generate  electric  polarization  under  an  inhomogeneous  mechanical  field,  and  yield  mechanical 
strain  under  an  inhomogeneous  electric  field.  In  1964,  Kogan  first  discussed  electric  polarization 
in  a  symmetric  crystal  by  an  inhomogeneous  deformation,  and  introduced  the  concept  of 
flexoelectricity  [1].  Experimental  and  theoretical  studies  have  since  confirmed  its  presence  in 
both  crystalline  and  amorphous  dielectrics  [17]-[23].  Recently,  it  has  been  reported  that  the 
scaling  effect  of  direct  flexoelectricity  has  great  potential  to  be  an  alternative  to  piezoelectricity 
for  micro/nano  devices,  suggesting  the  flexoelectric  effect  could  be  more  prominent  at  a 
micro/nano  scale  than  the  piezoelectric  effect,  even  in  state-of-the-art  piezoelectric  materials  [2], 
[24],  [25]. 


Although  there  have  been  extensive  studies  on  flexoelectricity  under  stress  and  electric  field 
gradients  [3],  [4],  [10],  [26],  [27],  investigations  of  flexoelectricity  under  thermal  fields  have 
been  rather  limited.  Because  of  the  temperature-dependent  dielectric  permittivity  of  flexoelectric 
materials  [28],  thermal  effects  should  be  a  topic  of  further  research.  It  is  expected  that  a  thermal 
field  will  induce  inhomogeneous  dielectric  properties  inside  the  flexoelectric  materials  [29], 
which  will  induce  the  stress/strain  gradient  under  an  electric  field  through  a  converse 
flexoelectric  effect.  Tagantsev  investigated  a  thermal  polarization  effect  in  flexoelectricity  in  the 
1980s  in  which  a  fundamental  relationship  between  a  linear  response  of  the  polarization  to  a 
temperature  gradient  was  discussed  [30],  [31].  Tagantsev  explained  the  thermal  polarization 
effect  on  the  flexoelectricity  from  the  concept  of  a  quasi-equilibrium  free  energy  density. 
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In  this  report,  we  investigated  the  thermal  polarization  effect  where  the  temperature-dependent 
dielectric  properties  of  BST  were  given  by  the  Curie-Weiss  law.  Derivation  of  the  relationship 
between  the  thermal  field  and  the  flexoelectricity  was  accomplished  through  an  anal3^ical  model, 
which  was  experimentally  verified. 

Secondly,  based  on  the  existence  of  the  converse  flexoelectric  effect  in  materials,  BST 
micro/nano-cantilevers  with  dual  comb  electrodes  were  designed  and  simulated  using  finite 
element  analysis. 

Thirdly,  we  investigated  a  new  material,  BaTi03-Bi(Zni/2Tii/2)03  perovskite  ceramics  (BT- 
8BZT),  at  a  broad  range  of  temperatures.  A  fully  comparison  and  analysis  were  presented  in  this 
report. 

Furthermore,  a  metal/ferroelectrics/semiconductor  structure  was  recently  tested  for  its 
flexoelectricity.  The  unique  direction  mechanism  were  obtained  from  the  experiment. 

Last  but  not  the  least,  during  the  reporting  period  we  also  conducted  designing  new  smart 
materials  by  alternating  its  structure.  A  preliminary  model  was  established  for  the  circular  curved 
microplates,  towards  smart  meta-structures. 
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5.  Summary  of  the  results 


5.1  Thermal  gradient  induced  flexoelectric  effect  in  hulk  Bao.67Sro.33Ti03 

5.1.1  Theoretical  modeling  and  simulation 

Showcased  in  Fig.  1,  a  temperature  difference  along  the  thickness  direction  of  the 
material  will  result  in  an  electric  field  gradient  due  to  the  temperature-dependent  dielectric 
properties  of  BST.  This  indicates  that  strain  can  be  generated  following  the  converse 
flexoeleetric  effeet.  The  theoretical  analysis  of  flexoelectricity  in  a  thermal  field  with  a  linear 
distribution  ean  be  given  by  the  Curie-Weiss  law  [15].  The  dieleetric  eonstant  of  the  ceramics  at 
paraelectric  phase  is  written  as 


1  1  T-T, 

-  = - + - - 

£  £  ,  C 

peak 

r  =  ri  +  (r, 

where  e  is  the  relative  dielectrie  eonstant,  is  the  peak  dieleetrie  constant,  is  the  Curie 
temperature,  C  is  the  Curie-Weiss  constant,  and  Tj  are  the  bottom  and  top  temperature  of  BST, 
X  is  the  axis,  and  L  is  the  thiekness.  For  Bao.evSro.ssTiOs,  s^gak  =  18000,  T^=  21  °C,  and  C= 
1.2x10^  °C  [28],  [32],  [33]. 


(3) 

(4) 


High  Temp. (72) 


Low  Temp.(ri) 


Figure  1.  Temperature  dependent  dielectric  permittivity  and  electric  field  curves. 
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The  boundary  condition  requires  the  electric  displacement  (D)  to  be  consistent  along  the 
height  direction,  and  integration  of  the  electric  field  equals  the  total  applied  voltage  (V) 


^  \^peatf  / 

=  -^  +  ^[2(ri-r,)+Cr2-ri)}|LxD 


(5) 


dEi, 


From  the  above  equations,  the  average  electric  field  gradient  (■^)  can  then  be  written  as 


dE^  ^  E^-E^  _  V 


n  -  T, 


dx, 


L2  1 


^penk 


f^(2(7’i-rj  +  (r,-r0) 


(6) 


where  Ei,  Ei  are  the  electric  fields  at  the  bottom  and  top  of  BST,  respectively. 

A  finite  element  analysis  was  adopted  to  confirm  the  above  calculations  using  the 
COMSOL  Multiphysics  (COMSOL  Inc.,  Los  Angeles,  CA).  The  electric  field  and  electric  field 
gradient  distribution  are  plotted  in  Fig.  2.  The  dielectric  constant  was  set  to  be  an  expression  of 
the  temperature,  as  depicted  by  Eq.  (3).  A  IV  voltage  drop  was  applied  to  the  top  surface  while 
the  bottom  surface  was  set  to  be  ground.  A  steady  state  simulation  was  implemented  to  obtain 
the  stable  temperature  distribution  and  the  electric  field  distribution.  Taking  into  account  a  10 
mm  thickness  (L),  21  °C  as  the  bottom  surface  temperature,  and  101  °C  as  the  top  surface 
temperature,  an  amplitude  of  the  average  electric  field  gradient  of  1.714x  10"^  V/m^  can  be 
observed  from  the  results,  supporting  the  analytical  estimation  from  Eq.  (6). 
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Figure  2.  The  electric  field  (a)  and  electric  field  gradient  (b)  distribution  under  1  V  with  the 
temperature  of  21  °C  at  the  bottom  and  101  °C  at  the  top  surface. 

Given  a  constant  temperature  difference  across  the  sample  surfaces,  Eq.  (5)  yielded  that 
the  average  electric  field  gradient  is  conversely  proportional  with  the  square  of  the  thickness  L. 
This  unique  property  implies  the  great  enhancement  of  the  electric  field  gradient  as  the  structure 
size  shrinks  down.  Moreover,  a  great  difference  in  temperature  leads  to  large  electric  field 
gradient,  which  is  desired  for  a  converse  flexoelectric  effect.  The  converse  flexoelectric 
coefficient  of  BST  is  characterized  as/iiii=6  x  10'^®  mW[9],  [34].  The  converse  flexoelectric 
effect  can  be  employed  here  for  actuation  applications  with  Eq.  (2),  where  the  strain  can  be 
readily  amplified  with  a  large  electric  field  gradient. 


5.1.2  Experimental  setup 

The  experimental  setup  to  measure  the  thermal  gradient-induced  flexoelectric  response  is 
schematically  shown  in  Fig.  3.  The  setup  was  arranged  on  a  floating  optical  table  (Newport,  ATS, 
Irvine,  CA)  to  eliminate  vibrational  noise.  The  AC  voltage  was  generated  by  a  power  amplifier 
(Trek,  2220,  Lockport,  NY)  under  excitation  from  a  function  generator  (Tectronix,  AFG3101, 
Lake  Mary,  FL).  To  measure  the  small  flexoelectric  displacement,  the  bottom  surface  of  the 
sample  was  clamped  on  the  A1  rod  in  which  a  conductive  silver  epoxy  was  adopted  as  an 
adhesive  and  the  rigidity  of  connection  was  carefully  checked  out  after  curing  time  of  silver 
epoxy.  The  axial  deformation  was  measured  using  a  high  resolution  (<10  pm)  laser  vibrometer 
(Polytec,  OFV-5000,  Irvine,  CA)  and  a  lock-in  amplifier  (Stanford  Research  System,  SR830, 
Sunnyvale,  CA),  so  the  measured  displacement  was  phase-locked  with  the  driving  voltage.  In 
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other  words,  any  other  background  vibration  (e.g.  building  vibration  or  hot  pate  vibration  if 
existed)  were  filtered  out  and  only  the  flexoelectric  vibration  was  recorded.  The  applied 
sinusoidal  voltage  ranged  from  200  V  to  1000  V  with  a  frequency  of  5  Hz.  Heat  was  applied  to 
BST  samples  through  the  heat  conductor  (A1  rod)  of  a  hot  plate. 


Hot  plate  (=  heat  source) 


Function  generator 


Al  rod  (=  heat  conductor) :  high  thermal  conductKity  (250  W/mK) 
Radiant  barrier  (prevent  radiation/convection  effect  on  BST) 


6 

Tt 

Silver  epoxy 
(adhesion) 


Lock-in  amplifier 


BST 


Power  amplifier 


Electrode  (Au) 


Figure  3.  The  experimental  setup  to  measure  the  thermal  gradient  induced  flexoelectric  response. 


To  assure  the  equilibrium  thermodynamic  status  of  BST  samples,  the  temperature  of  the 
bottom  and  top  surfaces  of  the  sample  was  checked  every  5  minutes  until  there  was  no  further 
temperature  change  detected.  The  stabilized  temperature  difference  between  the  top  and  bottom 
of  BST  samples  was  measured  by  a  thermocouple  (Fluke,  179,  Everett,  WA)  and  validated  with 
the  ID,  steady-state  heat  conduction  model  [35], 


Tjx-)  -  T, 
T,-T„ 


/  h  \ 

cosh{'m(L  —  :t)}  +  sinh{'m(L  —  :t)} 

cosh(7nL)  + 


m 


hP 

kA^ 


(7) 


where  r„,  are  the  ambient  and  heat  source  temperature,  and  L,  h,  P,  k,  Ac  are  the  length, 
convective  heat  transfer  coefficient,  plate  perimeter,  thermal  conductivity,  and  cross-section  area 
of  BST  samples,  respectively.  Fig.  4  represents  the  schematic  of  ID,  steady-state  heat  conduction 
model  and  the  analytical  modeling  result  under  the  condition  (7b=  100  °C,  Too  =  25  °C) 
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Figure  4.  Schematic  of  ID,  steady-state  heat  conduction  model  (left),  analytical  graph  (right). 

The  electrostrictive  effect  may  contribute  to  the  total  displacement,  in  which  the 
electrostrictive  displacement  is  featured  with  a  higher  frequency  signal  (2co)  than  the  exciting 
signal  (co)  [26].  Thus,  electrostrictive  strain  can  be  ruled  out  by  choosing  the  reference  frequency 
(oo)  of  the  lock-in  amplifier  to  be  the  same  as  the  exciting  signal  (oo).  Concerning  piezoelectricity 
in  BST,  when  the  BST  sample  temperature  is  above  the  Curie  temperature  (21  °C),  the  material 
should  theoretically  not  exhibit  any  ferroelectricity.  Nevertheless,  when  the  temperature  is  close 
to  the  Curie  temperature,  a  weak  persistence  of  the  macro-ferroelectric  region  may  exist  due  to 
local  nano-domains  [36].  In  Fig.  4,  a  small  displacement  (~2  pm)  was  observed  even  in  the  non- 
thermal  gradient  condition  (room  temperature,  23  °C),  which  reflected  the  residual 
ferroelectricity.  However,  compared  with  the  displacement  from  the  thermal  gradient-induced 
flexoelectricity,  the  piezoelectric  contribution  to  the  total  displacement  can  be  negligible.  As  the 
temperature  increases  (>  40  °C),  this  contribution  is  supposed  to  be  further  suppressed. 

BST  materials  used  in  this  work  are  ceramics  fabricated  by  Texas  Instrument  as  mentioned  in 
Ma’s  paper^^.  The  fabrication  process  involves  a  complex  wet  chemical  route  to  ensure  the  local 
and  global  homogeneity.  Four  BST  ceramic  samples  were  prepared  using  dicing  and  lapping 
processes,  and  the  dimensions  are  given  in  Table  I.  It  is  important  to  note  that  all  the  samples 
were  cut  from  the  same  BST  plate  along  one  direction  in  order  to  ensure  the  homogeneity  of  the 
samples.  Using  the  e-beam  evaporation,  Au  was  deposited  on  the  top  and  bottom  surfaces  of  the 
sample  to  a  200  nm  thickness. 
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Table  1.  Dimension  of  BST  samples. 


Sample 


Cross-section  Dimension  (ITi  x  IT2  x  L  mm^) 
shape 


BST  I 


Rectangle 


3  X  0.8  X  3 


BST  II 


Rectangle 


2  X  0.8  X  2 


BST  III 


Rectangle 


2  X  0.8  X  0.2 


BST  IV 


Rectangle 


1  X  0.8  X  0.1 


5.1.3  Results  and  discussion 

The  measured  face  displacements  of  samples  BST-I  and  BST-II  as  a  function  of  the 
applied  voltage  and  temperature  are  shown  in  Figs.  4(a)-(b),  respectively.  In  the  tests,  the  applied 
voltages  ranged  from  200  V  to  1000  V  at  5  Hz,  and  temperatures  at  the  bottom  of  BST  ranged 
from  40  °C  to  80  °C.  Displacement  of  BST  increases  directly  with  the  applied  voltage,  and 
shows  a  linear  relationship,  which  is  consistent  with  the  analytical  expectation.  Furthermore, 
compared  with  Figs.  4(a)  and  4(b),  under  the  same  conditions  of  applied  voltage  and  temperature 
range,  the  BST  sample  with  a  short  length  (2  mm)  shows  a  higher  displacement  than  that  of  the 
larger  BST  sample  (3  mm),  which  also  well  explains  the  feasibility  of  the  analytical  model.  It  is 
worth  noting  that  displacement  of  BST  is  not  proportional  to  the  temperature  difference  between 
the  top  and  bottom  of  the  samples,  even  though  the  higher  temperature  difference  clearly 
contributes  to  the  increase  of  electric  field  gradient  as  stated  in  Eq.  (6).  This  can  be  elaborated  by 
following  Eq.  (8), 


(8) 


Based  on  Eq.  (8),  the  total  displacement  of  samples  was  calculated  by  considering  the  converse 
flexoelectric  coefficient  and  the  electric  field  gradient.  The  converse  flexoelectric  coefficient, 
/ijtei’  shows  the  temperature-dependent  property.  Specifically,  decreases  as  the  temperature 

of  BST  increases  above  the  Curie  temperature  and  finally  becomes  zero  at  a  certain  temperature 
(~  101  °C)  as  shown  in  Fig.  6  [15], [28], [37].  Therefore,  despite  the  large  difference  in 
temperature  between  the  ends  of  the  BST  samples,  caused  by  the  high  temperature  applied  at  the 
bottom  of  the  BST,  further  increases  the  electric  field  gradient,  the  converse  flexoelectric 
coefficient  decreases  at  such  elevated  temperatures,  resulting  in  a  decrease  in  the  total 
displacement  of  the  samples  in  the  end  due  to  overweight  of  flexoelectric  coefficient  contribution. 
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Figure  5.  Displacement  as  a  function  of  applying  voltages  and  temperature  gradients 

(a)  L=2  mm,  (b)  L=3  mm 


Figure  6.  The  temperature-dependent  property  of  the  converse  flexoelectric  coefficient. 

To  further  investigate  the  scaling  effect  of  flexoelectricity,  all  samples  were  tested  at  400 
V  and  40  °C  bottom  temperature.  Fig.  7  represents  the  calculated  (based  on  Eq.  (2)  and  (6))  and 
measured  (based  on  the  measured  displacements)  strains  with  different  BST  samples,  in  which 
the  smallest  sample  (BST-IV)  shows  the  highest  strain  value  even  in  the  lowest  temperature 
gradient  condition.  It  also  shows  that  the  theoretical  analysis  is  well  consistent  with  the 
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experimental  results  on  a  basis  of  the  trend  of  data  even  though  there  are  some  discrepancies 
between  them. 


L  (mm) 


Figure  7.  The  calculated  and  measured  strains  of  BST  samples  having  different  lengths  and 

temperature  gradients. 
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5.2  Design  of  flexoelectric  nano-cantilevers  with  dual  comb  electrodes 
5.2.1  Electric  field  gradient  generation 

In  macro-scaled  devices,  the  electrical  potential  of  parallel  electrodes  only  would  generate  a 
perpendicular  electrical  field  distribution.  However,  in  micro-/nano-  scale,  according  to  the 
fringe  effect,  the  electrical  potential  applied  (+3  V)  on  one  unit  cell  structure,  in  Fig.  8,  would 
generate  electrical  fields  both  parallel  (x  direction)  and  perpendicular  (z  direction)  to  the 
electrode  at  the  edge  of  nanopatterned  electrodes.  Due  to  the  non-uniform  distribution  of  these 
electrical  field  in  the  space  (within  the  material),  the  electric  field  gradient  are  generated  in  both 
direction. 


Time=0.1  s  Surface:  Electric  potential  (V) 


Length  (nm) 

Figure  8.  Schematic  of  dual  comb  electrodes  beam.  (Top)  Applied  electric  potential  in  one  unit 

cell  of  dual  comb  electrodes. (Bottom) 

The  simulation  result  (COMSOL)  in  Fig.  9  confirmed  the  generation  of  electrical  field  gradient 
in  X  direction. 
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Time=0.1  s  Surface:  d(es.Ex,x)  (V/m^) 


Figure  9.  The  electrical  field  gradient  distribution  in  x  direction. 

According  to  the  conserved  flexoelectric  effect,  a  deformation,  denoted  as  S  in  Eq.  (9),  would  be 
generated  by  the  electrical  field  gradient. 

_  Ml  I  Ai2 

Cjj  dx  Cjj  dz 
//jj  =115[//C/  m] 

//,2  =  100[//C  /  m] 
c,i  =22.67  xlO“’[A/m"] 

//jj  and  //^2  are  the  flexoelectric  constant,  respectively,  is  the  elastic  stiffness  constant.  In  the 
case  of  study,  BST  is  used  in  further  discussion. 

5.2.2  Sample  design  illustration 

The  schematic  of  the  sample  design  is  illustrated  in  Fig.  10.  Both  the  top  and  the  bottom 
electrodes  of  the  beam  are  comb  electrodes.  The  beam  is  composed  by  two  separate  buck  BST 
samples.  The  insulated  ground  is  at  the  center  layer  of  the  beam,  which  would  provide  a 
protection  of  the  device  while  applying  the  electric  voltage.  Geometrical  parameters  of  the 
electrodes  and  the  beam  are  labeled  in  Fig.  10  as  well.  Based  on  the  result,  other  parameters  of 
thickness,  electrode  width,  and  electrode  gap  width  have  been  discussed  in  the  simulation  in 
order  to  look  for  an  optimized  parameter  set  in  further  fabrication. 
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Figure  10.  The  schematic  side  view  of  the  dual  comb  electrodes  beam  (4  unit  cells). 

5.2.3  Deformation  simulation 

A  finite  element  analysis  was  adopted  to  confirm  the  electric  field  gradients  using  COMSOL 
Multiphysics  (COMSOL  Inc.,  Los  Angeles,  CA).  The  electric  field  distribution  is  plotted  in  Fig. 
11.  And  the  electric  field  gradient  distribution  on  the  top  surface  and  the  bottom  surface  are 
plotted  in  Fig.  12.  The  boundary  condition  was  +3  V  applied  voltages  on  the  top  and  the  bottom 
electrodes.  The  simulation  results  show  that  the  electrical  field  and  field  gradient  are  non- 
uniform  in  the  sample,  which  agreed  with  the  previous  section. 


Time=0  s 


Length  (nm) 


Figure  1 1 .  The  electrical  field  distribution  in  x  direction. 
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Length  (nm) 


Figure  12.  The  electrical  field  gredients  distribution  in  x  direction. 

Furthermore,  based  on  the  simulated  value  of  the  electrical  field  gradient,  the  deformation  of 
BST  can  be  calculated  with  Eq.  (9).  We  evaluated  7  different  sets  of  parameters  to  discuss  the 
trend  of  increasing  the  deformation.  The  top  layer  contraction  and  the  bottom  layer  elongation 
are  summarized  in  Table  2.  Then  by  changing  the  electrode  width,  the  optimal  dimension  of  the 
dual  comb  electrodes  can  be  determined,  which  would  give  the  largest  top  contraction  and 
bottom  elongation  value,  listed  in  Table.  2.  Set  1-3  were  intended  to  evaluate  the  thickness.  As 
the  result,  thinner  ceramic  results  in  larger  deformation,  which  matched  the  scaling  effect  of 
flexoelectric  cantilevers  [8].  Set  3-5  were  intended  to  evaluate  the  gap  width  with  constant 
electrode  width.  It  turned  out  that  the  ideal  ratio  of  the  electrode  width  to  the  gap  width  would  be 
1.  The  fringe  effect  may  vanish  when  the  two  electrodes  get  closer.  Set  4,  6  and  7  were 
evaluating  the  electrode  width  with  the  constant  gap  width.  Narrower  the  electrode  width  is, 
larger  the  deformation  gets.  However,  in  further  experimental  study,  the  limitation  of  fabrication 
would  be  critical  to  select  the  desired  parameter.  These  evaluation  would  assist  the  selection. 


Table  2.  Calculated  deformation  based  on  different  ceramic  thickness,  electrode  width  and  gap 

width  (in  one  unit  cell) 


# 

Set 

Ceramic 
thickness  [nm] 

Electrode  width 

[nm] 

Gap  width 
[nm] 

Top  contraction 
[E-4  a.u.] 

Bottom  elongation 
[E-4  a.u.] 

1 

1600 

300 

600 

2.31 

2.28 

2 

800 

300 

600 

6.59 

6.57 

3 

400 

300 

600 

16.75 

16.72 

18 


4 

400 

300 

300 

19.93 

19.90 

5 

400 

300 

150 

18.55 

18.42 

6 

400 

30 

300 

45.79 

45.67 

7 

400 

150 

300 

28.24 

28.15 

5.3  Thermal  dependence  of  flexoelectricity  in  BT-8BZT 
5.3.1  Experimental  processes  and  results 

BT-8BZT  ceramics  is  of  interest  because  of  its  high  permittivity  at  a  broad  temperature  range. 
However,  the  flexoelectricity  of  BT-8BZT  has  never  been  examined.  In  this  study,  flexoelectric 
cantilever  beams  were  fabricated  with  BT-8BZT  and  BST  ceramics.  The  dielectric  permittivity 
was  examined  for  BT-8BZT.  Then  at  room  temperature,  the  flexoelectricity  of  BT-8BZT  was 
measured  by  using  a  bending  method,  while  piezoelectricity  was  directly  measured.  Furthermore, 
at  different  temperatures,  the  flexoelectric  transverse  coefficients  of  BT-8BZT  were  determined 
and  compared  with  BST  ceramics. 

Bulk  polycrystalline  samples  of  BT-8BZT  ceramics  were  prepared  using  a  conventional  solid 
state  synthesis  technique  by  Cann’s  research  group  at  Oregon  State  University,  as  described  in 
detail  by  Triamnak  et  al.*^.  The  sintered  pellets  were  cut  into  ceramic  bars  with  an  approximate 
size  of  1x5x10  mm^  using  a  diamond  wire  saw  (Well  Diamond  Wire  Saws,  Inc.,  Norcross,  GA). 
The  cantilever  beams  of  BT-8BZT  were  prepared  using  a  dicing  process.  Parameters  of  the  beam 
are  presented  in  Table  3.  Both  bottom  and  top  electrodes  consist  of  100  nm  (Au)/5  nm  (Ti) 
prepared  by  the  e-beam  evaporation  (Kurt  Lesker  Electron  Beam  Evaporation  System,  Jefferson 
Hills,  PA). 


Table  3.  Geometric  parameters  of  BT-8BZT  cantilever  beam. 


Thickness 

Length 

Width 

780  pm 

6.5  mm 

4  mm 

The  cantilever  set-up  have  been  discussed  in  the  previous  work  done  by  our  group  [2]  .In  addition 
to  the  room  temperature  set-up,  an  extra  thermal  chamber  built  with  the  ceramic  fiber  blanket  is 
inserted  in  Fig.  13.  To  record  the  temperature  inside  the  chamber,  both  the  thermal  couple  (TC) 
and  the  infrared  thermal  meter  (TM)  have  been  used.  After  3  minutes  of  blowing  hot  air  into  the 
chamber  by  a  heat  gun,  the  temperature  and  the  dropping  time  have  been  recorded  in  Fig.  14. 
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Figure  13.  Schematic  view  of  the  experimental  set-up. 


Figure  14.  Plot  of  the  temperature  (TC  &  TM)  vs.  the  cooling  time. 

The  thermometer  measures  the  temperature  through  the  infrared  radiation,  which  can  be 
absorbed  from  everywhere  in  the  insulated  environment.  This  would  lead  to  an  accumulation  in 
temperature,  which  would  be  higher  than  the  one  directly  measured  by  thermocouple.  Hence,  to 
simplify  the  set-up,  in  the  further  test,  only  the  thermocouple  is  used  to  record  the  temperature. 

Both  BST  and  BT-8BZT  cantilever  beam  had  been  tested  from  room  temperature  to  200  °C.  The 
flexoelectric  coefficient  is  calculated  in  Eq.  (10).  I  is  the  current,  f  is  the  frequency.  5  is  the  tip 
displacement.  A  is  the  electrode  area,  Lo  is  the  cantilever  length,  x  is  the  distance  between  the 
electrode  and  the  clamped  end. 


iKf  -A  35{I^-x) 


(10) 
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The  //j2  of  BST  and  BT-8BZT  are  plotted  in  Fig.  15.  Both  BST  and  BT-8BZT  have  high 

flexoelectric  coefficients,  25  pC/m  and  22  pC/m  at  25  °C,  respectively.  In  BST,  an  increase  in 
temperature  brings  the  fiexoelectricity  down  rapidly,  while  the  flexoelectric  coefficient  of  BT- 

8BZT  changes  more  gradually.  As  the  result,  the  flexoelectric  coefficient  of  BT-8BZT 

remains  around  12  pC/m  at  200  °C,  while  of  BST  is  less  than  3  pC/m  at  200  °C.  Since  the 
temperature  range  is  above  the  Curie  temperature,  weak  macro-ferroelectric  regions  may  exist 
due  to  the  presence  of  local  nano-domains  which  gradually  disappear  as  the  temperature  is 
increased.  To  verify  the  disappearance  of  this  residual  ferroelectricity,  the  direct  piezoelectric 
response  of  BT-8BZT  was  measured  over  the  same  temperature  range  and  the  data  is  shown  in 

Fig.  16.  As  indicated  in  the  figure,  the  piezoelectric  coefficient  ^^3  decayed  rapidly,  which 
indicates  that  the  piezoelectric  effect  eventually  disappears  in  the  paraelectric  phase.  Thus,  the 
flexoelectric  contribution  at  high  temperatures  dominates  the  mechanical-electrical  coupling. 


20  40  60  80  100  120  140  160  180  200  220 

Temperature  (°C) 


Figure  15.  Plot  of  //j2  of  BST  and  BT-8BZT  as  a  function  of  temperature. 
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Temperature  [°C] 


Figure  16.  Temperature  dependence  of  the  piezoelectric  coefficient 


^33  ofBT-8BZT. 


5.3.2  Discussion 

The  gradual  decline  of  flexoelectric  coefficient  with  increasing  temperature  in  BT-8BZT 
was  showed  in  this  work,  confirming  the  advantage  of  a  high  dielectric  permittivity  of  BT-8BZT 

in  paraelectric  phase  [38].  The  measured  flexoelectric  coefficient  of  BT-8BZT  is  25  pC/m  at 
25  °C  and  can  remain  above  12  pC/m  at  temperatures  up  to  200  °C,  which  suggests  potential 
applications  for  high  temperature  micro/nano-sensing.  Other  compositions  in  the  BiMeOs- 
BaTiOs  system  exhibit  temperature  stable  permittivity  over  a  much  wider  range. These  findings 
suggest  that  these  materials  may  be  more  effective  as  flexoelectric-based  sensors  over  a  broader 
temperature  range. 

5.4  Flexoelectricity  in  Semiconductor 
5.4.1  Theoretical  analysis 

From  the  beginning  of  the  study,  dielectric  insulators  have  been  the  objects  of  flexoelectricity. 
Recently,  semiconductors  have  been  studied  for  electrical  response  under  a  strain  gradient. 
Catalan’s  group  in  Spain  has  demonstrated  by  doping  single  crystals  of  wide-bandgap  oxides  to 
increase  their  conductivity,  their  effective  flexoelectric  coefficient  was  increased  by  orders  of 
magnitude.  Its  value  is  close  to  1  mC/m  [39].  The  samples  of  BTO-6  are  prepared  by  symmetric 
annealing.  The  oxygen  depletion  process  has  formed  a  Schottky  barrier  at  each  interfacial  barrier 
layer. 
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The  material  we  used  in  this  work  is  Pt/BTO/Nb:STO,  where  Nb:STO  is  also  an  n-type 
(vacancy)  semiconductor.  The  most  fundamental  characteristic  of  rectifying  metal- 
semiconductor  junctions,  such  as  Nb:STO,  is  the  Schottky  barrier  height  (SBH),  which  is  the 
energy  discontinuity  between  the  Fermi  level  of  the  metal  and  the  conduction  band  minimum 
(valence  band  maximum)  of  the  n-type  (p-type)  semiconductor.  Structures  like  Pt/BTO/Nb:STO 
is  classified  as  metal/ferroelectric  semiconductor  ferroelectric  tunnel  junctions(FTJs).  The 
resistive  switching  mechanism  appears  in  such  structure  with  different  polarization  directions.  In 
our  work,  the  barrier  direction  would  be  concerned  since  it  will  affect  the  strain  gradient  induced 
electrons  flow.  The  electrical  response  from  bending  in  different  directions  will  be  different  due 
to  the  different  resistivity. 

5.4.2  Experimental  results 

The  single -crystalline  0.7wt%  Nb-doped  SrTiOs  (Nb:STO)  is  chosen  as  the  semiconductor 
substrate.  A  100  nm  thick  BTO  film  was  deposited  epitaxially  on  (OOl)-oriented  Nb:STO 
substrates  by  pulsed  laser  deposition  in  a  typical  layer-by-layer  growth  mode.  Pt  top  electrodes 
were  deposited  on  BTO/Nb:STO  heterostructures  to  form  FTJs.  The  electrode  output  were 
monitored  from  the  Pt  electrode  and  the  Nb:STO  substrate  was  always  grounded  through  an 
indium  ohmic  contact  pad. 

Similar  experiment  set-up  was  used  as  described  in  5.3,  with  a  temperature  control  unit.  The 
sample  is  clamped  for  one  test  with  Pt  electrode  being  faced  up  and  the  other  test  with  the 
Nb:STO  face  up,  where  we  aimed  to  testing  how  the  generated  currents  are  different  when  the 
strain  gradient  is  in  the  same  direction. 

Fig.  17  shows  that  the  polarization  induced  by  the  flexoelectric  effect  is  linearly  proportional  to 

the  strain  gradient,  where  the  slope  indicates  the  average  xj^e  effective  value  for 
forward  direction  (Pt  face  up)  is  calculated  to  be  616.93+74.26  pCIm  at  room  temperature,  while 

the  effective  is  473.45+34.33  //C/m  for  backward  direction  (bottom  up). 
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Strain  Gradient  [m'^] 

Figure  17.  Flexoelectric  polarization  as  a  function  of  the  applied  strain  gradient  at 
room  temperature.  (Pink  &  red:  Pt  face  up;  green  &  lime:  bottom  up) 

5.4.3  Discussion 

In  this  work,  the  flexoelectric  switch  mechanism  is  found  in  Pt/BTO/Nb:STO,  which  is  due  to 
the  resistive  switch  mechanism.  The  experimental  results  can  support  the  expectation  from  the 
different  height  at  two  interfaces  of  Schottky  barriers. 

In  should  be  noted  that  the  magnitude  of  the  effect  is  large  comparing  with  conventional 
flexoelectricity  (~  100  //C/m).  Here  we  only  examined  with  the  one  type  of  semiconductor, 
transition  metal  oxides  with  the  perovskite  structure.  The  depletion  space  charge  region  on  the 
Nb:STO  surface  is  unique.  To  draw  a  further  and  more  thorough  conclusion  for  FTJs  structure 
type,  other  materials  should  be  studied  in  the  future. 

5.5  Flexoelectric  metamaterials 

5.5.1  Theoretical  analysis 

The  application  of  flexoelectricity  depends  on  the  structure  design.  The  circular  microplate 
structure  has  become  popular  in  the  design  of  smart  material  components.  Li  [40]  analyzed  the 
strain  gradient  elasticity  effect  on  a  flexoelectric  circular  microplate.  Furthermore,  an  enhanced 
flexoelectricity  in  curved  PZT  circular  plate  was  reported  by  Chu’s  group  [41,  42].  The  enhanced 
magnitude  of  the  mechanical-electrical  coupling  was  claimed  due  to  the  inconsistent  strain 
distribution,  which  is  caused  by  conversion  of  chemical  components  during  the  chemical 
reduction.  In  this  work,  an  analytical  model  was  studied  to  theoretically  investigate  the 
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flexoelectric  responses  in  the  curved  circular  plate  in  order  to  introduce  such  materials  to  sensing 
and  actuating  structures. 

Based  on  the  flexoelectric  theory,  a  flexoelectric  curved  circular  microplate  model  is  established. 
The  governing  equations,  boundary  conditions  and  initial  conditions  are  derived  following 
Hamilton’s  principle. 


Figure  18.  Schematic  of  a  flexoelectric  curved  circular  microplate  and  its  coordinate  system 

As  shown  in  Fig.  17,  the  thickness  of  the  curved  circular  microplate  is  h.  The  curvature  radius  of 
the  central  plane  is  R.  The  central  angle  is  .  The  curved  circular  microplate  is  subjected  to  a 
symmetric  distributed  force  q  and  a  voltage  V  on  the  upper  surface.  A  spherical  coordinate 
system  is  used  in  the  microplate  with  the  6-(p  plane  coincident  with  the  central  plane  and  the  r- 
axis  along  the  curvature  radius  direction.  According  to  the  Kirchhoff-Love  thin  shell  theory,  the 
displacement  components  of  the  curved  circular  microplate  in  the  spherical  coordinate  system 
are  taken  as 


Uy.=w{(p,t),  u„  = 


R 


■  +  1 


Uff  =0, 


(11) 


Here,  z  =  r-R  and  u^,  u^,  are  the  displacements  in  the  curvature  radius,  meridian  and 
latitude  direction,  respectively.  w[<p,t),  u{(p,t)  are  the  deflection  and  the  meridian  displacement 
of  central  plane,  respectively. 

Then,  we  get  the  flexoelectric -polarization  in  the  direct  flexoelectric  effect  as 
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Here,  rj^r^'neer^'nrzixp^'nree  denote  the  strain  gradients  components.  Ai22’A212  ^^e  the  flexoelectric 
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coefficients,  is  the  electric  potential  and  V^= — ^  +  cot(^)— . 

d(p  d(p 

In  the  simulation  of  the  direct  flexoelectric  effect,  the  induced  surface  charge  linear  density  is 


2(^£>)  =  2;ri?sin(^)C^^^ 


\2j 


(13) 


is  the  capacity  per  unit  surface  area  defined  as  /  h  and  is  the  dielectric  constant. 

And  hence  the  total  amount  of  collected  charges  is 


Qcou  =  ^m{(p)d(p 
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For  the  converse  flexoelectric  effect,  the  larger  deflection  is  usually  the  goal.  In  this  case,  the 
applied  voltage  V  is  found  to  be  approximately  equivalent  to  an  applied  uniform  pressure 


^(Ai122  '*'^^1212)^ 


combined  with  an  applied  boundary  moment  whose  linear  density  is 


A 122^ . 


To  theoretically  investigate  the  flexoelectric  responses  in  the  curved  circular  plate,  some 
numerical  simulations  and  discussions  are  given.  Here,  the  BST  ceramics  material  is  used,  whose 
material  parameters  are  listed  in  [4].  In  the  direct  flexoelectric  effect,  =;t/10  and  the  applied 
pressure  is  O.lMPa.  In  the  converse  flexoelectric  effect,  h  =  l//m  and  the  applied  voltage  is  20V. 
In  both  cases,  R  =  20h  /  and  the  clamped  boundary  is  considered. 


In  the  direct  flexoelectric  effect,  the  relation  between  the  bending-induced  charge  density  and  the 
thickness  is  shown  in  Figure  19(a).  In  the  converse  flexoelectric  effect,  the  relation  between  the 
maximum  voltage-induced  deflection  (deflection  at  the  center)  and  the  center  angle  is  shown  in 
Figure  19(b). 
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5.5.2  Analytical  results 


Figure  19.  (a)  Direct  flexoelectric  effect;  (b)  converse  flexoelectric  effect 

As  shown  in  the  Figure  19,  in  the  direct  flexoelectric  effect,  the  bending-induced  charge-density 
(A  is  the  surface  area)  is  found  significantly  increasing  as  the  thickness  decreases  to  nanometer 
scale,  which  demonstrates  the  size-effect  of  the  flexoelectricity  in  curved  circular  plate.  In  the 
converse  flexoelectric  response,  voltage-induced  bending  will  be  generated  in  the  curved  plate 
rather  than  in  the  flat  circular  plate,  which  demonstrates  some  advantages  for  curved  circular 
plate  over  the  flat  plate.  Here,  when  the  central  angle  tends  to  zero,  the  curved  plate  gradually 
turns  to  be  a  flat  circular  plate. 

5.5.3  Discussion 

The  theoretical  model  of  a  flexoelectric  circular  plate  is  established  to  calculate  its  direct  and 
converse  flexoelectricity.  Based  on  this  theory,  we  are  able  to  determine  the  structural  sizes 
according  to  the  target  responses.  However,  the  calculated  sizes  may  be  only  used  as  reference 
values  since  many  other  factors  in  the  experiment  such  as  fabrication  precision,  test  precision  etc. 
will  also  affect  the  results. 

5.6  Conclusion 

During  the  final  reporting  period,  the  following  key  results  were  obtained. 

(1)  The  flexoelectricity  of  Bao.eTSro.ssTiOs  (BST)  was  investigated  by  generating  temperature 
gradients  along  the  lengths  of  samples  with  symmetric  geometry.  An  electric  field 
gradient  induced  by  a  thermal  gradient  was  analyzed  based  on  the  temperature-dependent 
dielectric  property  of  BST.  The  strain  was  then  experimentally  verified  due  to  the  electric 
field  gradient.  Experimental  results  suggest  converse  flexoelectric  effect  of  BST  samples 
with  symmetric  geometry  in  a  thermal  field.  This  result  was  not  only  consistent  with  the 
theoretical  prediction,  but  it  also  followed  the  scaling  effect  of  flexoelectricity. 
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(2)  The  converse  flexoelectric  effect  was  investigated  in  the  BST  beam  structure.  The  dual 
comb  electrodes  has  induced  the  electrical  field  gradient  and  increase  the  beam  deflection. 
The  grounded  insulation  has  simplified  the  further  fabrication  work  by  eliminating  the 
alignment  between  the  top  and  bottom  electrodes. 

(3)  The  novel  material  BT-8BZT  has  been  studied  for  its  flexoelectricity.  The  smooth 
variation  of  flexoelectric  coefficient  has  been  observed  in  the  experiment  in  the 
temperature  range  from  25  °C  to  200  °C.  The  thermal  dependence  of  flexoelectricity  in 
BT-8BZT  can  be  further  adopted  in  high  temperature  sensing  research. 

(4)  The  FTJ  metal/ferroelectrics/semiconductor  structure  has  been  tested  for  its 
flexoelectricity.  The  resistive  switch  effect  not  only  has  a  significant  in  the 
semiconductor  field,  but  also  opens  a  new  direction  of  design  for  smart  materials  in 
sensing  and  actuating  applications. 

(5)  The  theoretical  model  for  flexoelectric  metamaterials  has  been  derived  towards  deeper 
understanding  and  verification  of  its  intrinsic  principle.  However,  up  until  now,  only  a 
single  layered  circular  curved  microplate  has  been  tested  under  the  analytical  model.  It  is 
still  too  early  to  draw  a  firm  conclusion  from  the  current  simulation  results  for  layered 
curved  microplates. 
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